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The general features of Microcanonical Thermodynamics (MT) as applied to the 
fragmentation of hot nuclei and atomic clusters are discussed. MT is the most 
fundamental form of any thermodynamics since Boltzmann. With modern compu- 
tational techniques it is for the first time possible to explore it in nontrivial cases. 
First and second order phase transitions in finite systems can unambiguously be 
identified by the caloric equation of state T(E/N). The three characteristics of 
phase transitions: The transition temperature, the latent heat and the interphase 
surface tension can be well defined and calculated in finite systems. Against com- 
mon believe and in contrast to conventional (canonical) thermodynamics, MT of 
relatively small systems reflects the thermodynamic behavior of bulk systems in 
great detail and surprising accuracy. By the Laplace transform from the micro- 
canonical to the canonical ensemble many of the signatures of phase transitions 
are smeared out or get lost. Also some of our believes about phase transitions 
must be corrected: Microcanonical phase transitions show distinct peculiarities 
very different from conventional ph-tr.: E.g. the specific heat is negative. Phase 
separation is usually unsuitable to identify ph-tr. in small systems and in sharp 
contrast to the caloric equation of state. In contrast to ordinary (canonical) ther- 
modynamics MT allows the system to become inhomogeneous or to fragment, at 
first order phase transitions it gives insight into the coexistence region. Here the 
form of the specific heat c{E/N) connects transitions of first and second order in 
a natural way. The "phase transition" towards fragmentation is introduced. The 
similarities and differences to the boiling of macrosystems are pointed out. 

Nuclear friction is likely responsible for the ergodic expansion of a multifragmented 
nuclear source up to freeze-out densities of ~ 1/6 normal density. Various recent 
experimental results which seem to contradict our scenario are even confirming it 
when investigated more thoroughly. The fragmentation "phase-transition" (multi- 
fragmentation) in nuclei is by many reasons not the same as the liquidi-cas transi- 
tion in nuclear matter. The main reasons were already published 1984^1 



1 Introduction 

Our definition of Microcanonical Thermodynamics (MT) uses only mechanics. 
The topology of the total accessible N-body phase space reflects (or implies?) 
the behavior of many interacting N-body systems because their dynamics is 

"We use quotation marks to emphasize the limited analogy to standard phase transitions, 
CRIS96, Catania May 27-31,1996. "Critical Phenomena and Collective Observables" 
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often chaotic. Then the dynamical evolution of many replica of the same 
system under identical macroscopic initial conditions follows the structure of 
the underlying N-body phase space. It is ergodic. In nuclear fragmentation 
this is presumably due to the strong and short ranged friction between moving 
nuclei in close proximity. Friction between atomic clusters is yet unknown but 
quite likely it exists there also. 

First we have to discuss the concept of thermodynamics of small systems 
in general and especially of phase transitions. The first question is from which 
size on do phase transitions exist in small systems? How can one define them? 
Is a cluster of ~ 100 particles big enough? We will show that this is possible 
and one can unambiguously distinguish continuous (second order transitions) 
from discontinuous (first order) transitions by the form of the caloric equation 
of state Tthd{E). Before proceeding further it is important to realize that 
isolated nuclei or clusters must be treated microcanonically. Usually there is no 
external heat- or particle bath which defines the temperature, the pressure, or 
the chemical potential. Microcanonical Thermodynamics is the proper theory 
for isolated small systems. 

A microcanonical ensemble has some peculiarities: It does not have a 
positive definite heat capacity. In fact at a phase transition of first order the 
specific heat c(e) becomes negative in general. Therefore the classical signal of 
a peak in the specific heat is not useful to characterize a phase transition in 
small systems. 

The second peculiarity of Microcanonical Thermodynamics is that it allows 
the system to become inhomogeneous: At first order phase transitions several 
regions of one phase coexist with other regions of the other phase. MT allows 
diiferently to conventional thermodynamics the coexistence of regions with 
different energy-density. The partitional entropy of the spatial fluctuations 
is an important part of the total entropy. Small many-body systems and also 
large systems under long-range forces have an important new structural "phase 
transition" which does not exist in infinite homogeneous systems: They may 
fragment into few relatively large pieces. Tycical examples are nuclear multi- 
fragmentation, as was predicted very earlyla'Ba, but also the fragmentation of 
atomic clusters. 

In this case the size fiuctuations of the fragments at the transition are of 
the order of the size of the system itself. Then one cannot ignore the "droplets" 
compared to the nucleonic or monatomic vapor anymore. The "phase tran- 
sition" is not determined alone by the equilibrium of the homogeneous liquid 
with the homogeneous gas, which in conventional (canonical) thermodynamics 
is controlled by the equality of the chemical potentials of liquid and gas as 
we are used to in conventional (grandcanonical thermodynamics). Often the 
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number of droplets is similar or larger than the number of free nucleons or 
monomers. This is one of the main lessons relevant for general physic which 
can be learned from nuclear fragmentations. 



2 Microcanonical Thermodynamics 

Microcanonical Thermodynamics explores the topology of the N-body phase 
space and determines how the volume fl^ of the accessible phase space — 
more precisely the number of quantum states — depends on the fundamental 
globally conserved quantities of total energy E — N * e, angular momentum 
L, mass (number of atoms) N, charge Z, linear momentum p, and last not 
least the available spatial volume V of the system. Thiajiefinition is the basic 
starting point of any thermodynamics since BoltzmannB 
The entropy is defined as the logarithm of 

S{E, V, N) = Ns{e = E/N) = ln{rL{E, V, N)) (1) 

and the thermodynamic temperature Ttud is defined by 

Tthd = (3) 

By Laplace transform of ^{E, V, N) one steps from the "extensive" variables 
like E, V, N to the intensive ones like T, P, (i. E.g. the Gibbs grand partition 
function and the Gibbs grand potential are then 

Z{(3,P,fi) = JJJ n{E,V,N)e-'^'^^+^^->'^^ dE dV dN, (4) 
G(AP,/i) = -Tln[Z{p,P,^^)]. (5) 

We take Boltzmann's constant k = 1. In the same way one gets the canonical 
partition function and the free energy as 

Z{I3,P,N) = JJ n{E,V,N)e-'^^^'+^''UE dV, (6) 
F(r,P,iV) = -Tln[Z{l3,P,N)]. (7) 

3 Differences between microcanonical and canonical ensemble 

According to van Hove a system of N particles interacting via short range 
two-body attractive forces with hard cores is thermodynamically stable, the 
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thermodynamic limit of N,V cxDj^^y^g exists for such svstems, intensive 
quantities hke the specific energy have finite limiting valuesEl Then, the ther- 
modynamics derived from the microcanonical partition sum Q{E, N,V) and 
the one derived from the canonical Z{P,P,N) or the grand canonical parti- 
tion function Z{P, P, ^) (usually) coincide. Outside of phase transitions of first 
order, the relative fluctuations /S.E/E, or Ag/g vanish cx 1/^/N. 

This is quite different at phase transitions (T = Ttr) of first order and for 
finite systems where the microcanonical and the (grand) canonical ensemble 
differ essentially. In the canonical ensemble the energy fluctuations Ae per 
particle remain finite even in the thermodynamic limit. ((Ae)^|Tt^ ^ Hat, the 
specific latent heat). Consequently, the difference between the microcanonical 
and the canonical ensemble persists at transitions of first order and we must 
expect both ensembles describe different physical situations. 

Systems interacting via long range forces like unscreened Coulomb or the 
centrifugal force when they are rapidly rotating don't have a thermodynamic 
limit and must be described by the microcanonical ensemble. Such systems 
fragment macroscopically into, in general several, regions of high density — 
condensed matter — and also into, in general several, regions of low density — 
vapor or may be empty space. Differently from conventional thermodynamics 
where systems which must be in a homogeneous phase at fixed temperature ev- 
erywhere, here the system is most likely inhomogeneous. The inhomogeneities 
and their fluctuations are more important in characterizing the state of the 
system than any mean values. In contrast to thermodynamics of the homoge- 
neous bulk, in small systems or large systems with long range forces the entropy 
connected to different partitions of the system is an important part of the total 
entropy. Familiar formulas like the one-particle entropy 



are useless for calculating the total entropy. 

4 First and second order phase transitions in small systems 

Macroscopic systems have a discontinuity in the specific heat Cbuik{T) at first 
order phase transitions. Cbuik{T) may have a finite peak at T « Ttr. On top 
of this there is a spike qiat5{T — Ttr). With finite resolution it shows up as 
jump in Cbuik{T) by the latent heat qiat- A typical example is the specific heat 
of bulk sodium, fig.(l) at the melting transition. In contrast, a transition of 
second order is continuous at the transition temperature where Cbuik{T) has 
(in the example of the Ising model) a logarithmic singularity in T — Ttr- 




(8) 
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Microcanonical Thermo- 
dynamics gives new and 
deep insight into this. It 
wiU further aUow to ex- 
tend the concept of phase 
transitions to systems not 
treated before by ther- 
modynamics hke systems 
with long rang forces or 
strongly rotating systems. 
We begin with the dis- 
cussion of microcanonical 
phase transitions in stan- 
dard model systems in 
which phase transitions of 
first and second order are 
well known. As example we take the 2-dimensional 10-states Potts anodel fpc 
which the asymptotic thermodynamics is even known analyticalljcl In rcfllj 
we determined the three basic parameters of phase transitions of first order 
within Microcanonical Thermodynamics, the transition temperature Ttr, the 
specific latent heat qiati and the specific interphase surface entropy l^Ssurf 
for a system with nearest neighbor couplings. It was demonstrated that for 
surprisingly small systems the values of these three parameters are closer to 
their asymptotic values than in the canonical ensemble. This is so because 
most of the finite-size scaling is due to the large, but trivial, exponent ip-itiije 
Laplace transform eq. from the micro- to the canonical partition sumEil'Ej 

The two types of phase transitions are distinguished by the form of the 
microcanonical caloric equation of state Ttudi^)'- A transition of first order 
has a backbending caloric equation of state T"^ — (3(e) c.f. figure (2b). For 
a system with nearest neighbor interactions the area between /3(e) and the 
"Maxwell" line (3 = 1/Ttr is twice the interphase surface entropy Ass„r/lij The 
left darkened area is the defect of entropy I^Ssurf — J^^ Pde that the system 
'pays' for introducing interphase surfaces , which it finally gets back when the 
whole system is converted to the new phase at e = £3, right darkened area, and 
the interphase surface disappearesEj In the bulk the transition is discontinuous 
as function of T or /3 and "jumps" from the liquid {s < Si) to the gas branch 
{e > £3) of the caloric curve, fig. (2). As a function of the specific energy e the 
transition is however continuous. With rising e the system passes smoothly 
from the liquid phase over a mixed phase with coexisting large fluctuations 
of the two phases ( "gas bubbles" and "liquid droplets" ) to the pure gas phase 
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when the specific energy is increased by the specific latent heat qiat- Figure (2a) 
shows the corresponding specific entropy s(e) = /q P{£')de' . The transition is 
characterized by the convex intruder in s{e) of depth Asgurf- Figure (2c) 
shows the specific heat capacity 

de -(3'^ 

as a function of the specific energy e. (Here numerical fluctuations in /3(e) in 
figure (2b) have been smoothed). One can see within the coexistence region of 
El < e < £3 (shaded area in fig (2) ), the microcanonical specific heat has two 
poles and becomes negative in between. 

As the convex intruder in the specific entropy s(e), fig. (2a), is forbidden 
by van Hove's theorem in the canonical ensemble for an infinite number of 
particles, conventional thermodynamics of the bulk is blind in this energy 
intervalB Here the Laplace transform eq.(TO-,has additional stationary points 
and the canonical bulk would be unstabletll It can only see the branches of 
c(e) in the regions e < ei and e > £3. Thus the canonical specific heat ci,uik{T) 
will be positive and approach finite values at the transition of first order. At 
T = Ttr ciiuik{T) has an additional peak = qiatS{T — Ttr)- 

If the specific latent heat qiat — > and the specific interphase surface 
entropy /S.Ssurf — *■ the caloric equation of state gets only a saddle point at 
the transition. Then E{T)/N as well as T{E/N = e) become single valued, 
the transition is continuous in the canonical as well as in the microcanonical 
ensemble. We have a phase transition of second order. The two poles of c{e) 
merge and c(e) or ci,uik{T) has a singularity at the transition point etr^Ttr- 
Both slopes of c(e) are fully accessible in the canonical treatment of the heat 
capacity. Consequently, from the caloric equation of state T{e) it is always 
possible to identify and distinguish both kinds of transitions. In Microcanonical 
Thermodynamics the relation between the two is very natural, transparent and 
simple. 

It is further instructive that in finite realizations of the two-dimensional 
Potts model with q = 10 spin orientations at each lattice point it was not 
possible to see a clean separation into a compact region of ordered spins and 
a compact region of disordered spins at energies inside the coexistence region 
even for a lattice of 100 * 100 points. There were always several "gas bubbles" 
and "droplets" fluctuating over the lattice and prohibiting large interphase 
surfaces. Nevertheless the caloric equation of state T{£) is already close to its 
asymptotic form. We can conclude from this observation that the other clas- 
sical signal of a transition of first order, a clear separation of the two phases, 
is not a useful signal of a transition of first order in small systems. 
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Fig. (2) Specific entropy s(e) — J^^ 0micro{^)de 
vs. the specific energy e for the 2-dim. Potts 
model with g — 10 on a 200 * 200 lattice. In or- 
der to visuahze the anomaly of the entropy the 
linear function a be (a — s{e — 0.17), b — 1.42) 
was subtracted. Because we use periodic bound- 
ary conditions one needs two cuts to separate the 
phases and the depth of the convex intruder is 
twice the surface-entropy. 

b) Inverse temperature ^microi^) — ^/T{e) as di- 
rectly calculated by MMMC 

c) Specific heat c{e) - -/3^ / {dl3/de). The canon- 
ical ensemble of the bulk jumps over the shaded 
region between the vertical lines at ei and £3. 
This is the region of the coexistence of two phases 
one with ordered spins, the other with disordered 
spins. Here c{e) has two poles and in between it 
becomes negative. The canonical thermodynam- 
ics is blind to this region. Observe that the poles 
are inside £1 < £ < £3, i.e. the canonical specific 
heat remains finite and positive as it should. 

5 Signals of a "phase transition" in nuclear fragmentation 

In the review article refEl I proposed the caloric equation of state T{e) for 
^^^Xe. This had two anomalies compared to the standard parabolic depen- 
dence of a Fermi-gas T cx ^/e*, fig. (13). One is quite pronounced with even a 
backbending of T(e*) at e ~ 2.5MeV, T ~ 4.5MeV and a second one is less 
pronounced and has no backbending is at £ ~ 5MeV, T ~ 6MeV. The first one 
had a width qiat ^ IMeV/nucleon and the second ~ 1.5MeV/nucleon and may 
even be a transition of second order{l). The latter one shall not be discussed 
here. This is likely the transition found recently at GSIolij In the above re- 
view article the first anomaly was compared to the apparent slope temperature 
Tapp of evaporated a-particles from "^^S+Ag and of ^^0+Ag from the Texas 
A&M group whic h show the same narrow anomaly at similar energies and 
temperaturestjO However, these data had too large error-bars to allow any 
firm conclusion. The interesting aspect of this first anomaly is that one may 
see it in theory as well as in experiment in a- evaporation spectra even though 
in the model they are linked to the « sudden opening of the IMF production. 
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The experiment by Chbihi et 
al. of incomplete fusion re- 
actions of 701 MeV 2^Si + 
^°°Mo gives fupLh^r evidence of 
this transitionliflllj We plot in 
fig. (3) £*, the excitation energy 
per nucleon vs. Tapp, where 
Tapp is the slope of the raw 
evaporation spectra for protons 
and alpha particles. The curves 
give the Tapp{s*) dependence 
deduced from the microcanon- 
ical statistical multifragmenta- 
tion model (MMMCV-Lisipg its 
standard parametersu3li3 Also 
the experimental uncertainties 
for the proton and alpha curves 
are given. The similarity of the 
shapes of the experimental and 
simulated Tapp{e*) for the a- 
spectra is quite evident. The 
differences between the shapes 
of these curves and the parabo- 
lic dependence (dotted curve) expected for a simple Fermi gas is clearly outside 
the experimental error margins indicating that some additional degrees of free- 
dom, which are apparently included in the [MMMC) model, become significant 
in this energy range. The proton data are not so clear. They do indicate a 
similar anomaly but are more close to the parabolic Fermi-gas form than to the 
MMMC curve. Maybe the protons diffuse too fast out of the expanding soup 
of fragments and do not explore the structure of the accessible phase space 
sensitively enough. 

The experimental analysis of the data provides the values of the mass A. 
charge Z, excitation energy E* , and angular momentum L of the sourceEj 
The freeze-out radius Rf was taken as its standard value of 2.2 A^/'^ fm, this 
means that we simulate a "phase transition" at constant volume. The results 
of MMMC calculations, performed by O.Schapiro, with these input values, 
were subjected to the same software filter as the experimental set-up which, 
most importantly, selects only those events with one big residue. The mass of 
the residue was chosen to be Ares > 90, which is close to Ares estimated from 




Tapp (McV) 

Fig. (3) Experimental and theoretical (with MMMC) 
caloric equation of state, Tapp(,e*) for p and a. The 
horizontal error-bars give the statistical unccrl^nty to 
extract the slope from the raw spectra in rcftJ. Dif- 
ferent methods to determine the excitation energy lead 
essentially to a parallel up or down shift of the curves 
by the amount indicated by the vertical bars at the low- 
est and highest data point. The dash-dotted curve is a 
Fermi-gas calculation (e = T^/8). 
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the experimental data (the mass of the residue could not be measured). 

The theoretical value of Tapp was extracted from fitting, as was done for 
the raw experimental spectra. Similar to the experiment the individual tem- 
peratures for the different particle species are slightly different in MMMC from 
the thermodynamic temperature Tthd cf. figure (3). This interesting detail is 
due to the intrinsic fluctuation of Tthd in the microcanonical ensemble — a 
result outside conventional thermodynamics. 

The values of i?/ and A^es do not influence the general shape of the 
7app(e*) curves. However, the Tapp{s*) curves shift along the Tapp - axis if 
different values of these parameters are used. The shifts produced by reason- 
able changes in Ares are larger than those produced by reasonable changes in 
Rf. We checked that the anomaly in T{e) is not due to the changes of the 
angular momentum from L — 18.2 to 48.8?i. It exists also at L = 0. 

While the similarity of the shapes of the experimental and simulated 
7app(e*) is quite reasonable for a, d, and t-particles, (see Fig. (3)) significant 
differences exist in detail. The simulated curves for d and t (not shown here) 
have the same S-shape as the data (and the a-s) but are shifted towards lower 
values of Tapp. The higher Tapp values of the experimental deuteron and triton 
spectra might again be an indication of a faster diffusion of these less bound 
fragments from a hotter stage of the early fragmented system. However, differ- 
ently to the protons their diffusion is on the other hand slow enough to realize 
the transition in the underlying structure of the pJiase space. This interpre- 
tation is supported by a recent experiment by refo showing a simultaneous 
early fragmentation and a smaller "decay" (diffusion) time for d and t. 

The experimental data also suggest an association between the onset of 
IMF production and the anomaly in rapp(e*). In both theory and experiment 
the multiplicity of Li-fragments just starts to rise at the same excitation. The 
success of the MMMC model in showing the same anomaly in the caloric 
equation of state in a nearly automatic manner has a potential far reaching 
consequence. The central assumption of this model is an ergodic mixing of 
the fragmented system filling the accessible N-fragment phase space with a 
well defined freeze-out region uniformly. This is quite different from sequential 
binary fragmentation, with independent motion of fragments outside the binary 
barrier. Here an anomaly in E{T) would not have such a natural explanation. 

Taking all evidence together: The anomaly in all four spectra (proton, 
deuteron, triton, and alpha) at the same excitation energy as predicted by 
MMMC and also the earlier data of the Texas A&M-group, in spite of their 
huge error bars there is in my opinion no doubt that this is a signal of a "phase 
transition" due to the relatively sudden opening of additional phase space at 
the onset of fragmentation, the lower one predicted in refE^. Due to the strong 
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stochastic (dissipative) couphng of the various fragmentation channels this 
"phase transition" is felt in all chajMuds even in the pure evapacation channels. 
The data of the Texas A&M grouji! the data of Chbihi et alEJ, as well as the 
results of the Aladin collaboratioicJ are the first^jxperimental caloric signals of 
a phase transition in a nucleus. The data of refllj have now smaller erxor bars 
and approach the higher (second order ?) transition mentioned abovell3 They 
show a slow monotonic increase of T with rising excitation. I.e. they do not 
anymore exhibit the dramatic plateau in T{E/A) with the spectacular "snecific 
latent heat" of ~ 5 MeV, close to the total binding energy per nucleouEj 

1. D.H.E. Gross. 1984 INS-RIKEN Symp.,J.Phys.Soc.Jap.54(1985)392 

2. D.H.E. Gross and Meng Ta-chung. In 4th Nordic Meeting on Intermedi- 
ate and High Energy Physics, page 29, Geilo Sportell, Norway, January 
1981. 

3. J. P. Bondorf. In C.H. Dasso, editor, Nuclear Physics, Proceedings of 
the Nuclear Physics Workshop, I.C.T.P., Trieste, Miramare, Italy, 5- 
30 October, 1981, pages 765-770, Amsterdam, New York, Oxford, 1982. 
North Holland. 

4. D.H.E. Gross, L. Satpathy, Meng Ta-chung, and M. Satpathy. 
Zeit. Physik, A 309:41, 1982. 

5. Ludwig Boltzmann. Vorlesung iiher Gastheorie. Number 1. Akademi- 
sche Druck-u. Verlagsanstalt, Graz, 1981. 

6. L. van Hove. Physica, 15:951, 1949. 

7. G. Borelius. Solid State Physics, 15:1, 1963. 

8. R. Hultgren, R.L. Orr, P.D. Anderson, and K.K. Kelley. Selected Values 
of Thermodynamic Properties of Metals and Alloys. John Wiley and 
Sons, New York, 1963. 

9. R.J. Baxter. /. Phys., C6:L445, 1973. 

10. D.H.E. Gross, A. Ecker, and X.Z. Zhang. Annalen der Physik, in print, 
1996. 

11. A. HiiUer. Z.Phys.B, 95:63-66, 1994. 

12. M. Promberger. preprint, Erlangen, 1996. 

13. D.H.E. Gross. Rep. Progr. Phys., 53:605-658, 1990. 

14. J. Pochodzalla et al. Phys. Rev. Lett., 75:1040-1043, 1995. 

15. C. Schwarz et al. GSI-Nachrichten, 01-96:5-16, 1996. 

16. G. Nebbia et al. Phys. Lett. B, 176:20-25, 1986. 

17. R. Wada et al. Phys.Rev.C, 39:497, 1989. 

18. A. Chbihi et al. Phys.Rev.C, 43:652;43:666, 1991. 

19. A. Chbihi, D.H.E. Gross, O. Schapiro, S. Salou, and L.G. Sobotka. sub- 
mitted to Phys. Rev. Lett., 1995. 

20. C.J. Gelderloos et al: Phys. Rev. Lett, 75:3082-3085, 1995. 



10 



